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THREE-DIMENSIONAL POTENTIAL FLOW OVER HILLS AND OVAL MOUNDS 

by Robert Siegel 
Lewis Research Center 


SUMMARY 

An analysis was made of the potential flow behavior for an initially uniform flow 
passing over a single axisymmetric hill, an oval mound, and a combination of two hills. 
Small perturbation theory was used, and the resulting Laplace equation for the perturba- 
tion velocity potential was solved by using either a product solution or a Green's function. 
The three-dimensional solution is of interest in calculating the pressure distribution 
around obstacles, the flow of pollutants carried by the wind, and the augmentation of wind 
velocity for windmill siting. The augmentation in velocity at the top of a hill was found 
to be proportional to the hill height relative to a characteristic width dimension of the 
hill. An axisymmetric hill produced about 20 percent less velocity increase than a two- 
dimensional ridge having the same cross-sectional profile. 


INTRODUCTION 

The behavior of three-dimensional flow over irregular terrain is of interest for de- 
termining the flow path of pollutants, the pressure distribution over obstacles, and the 
augmentation of the wind velocity for selecting windmill sites. Some examples of air 
pollution studies are given by references 1 and 2, while reference 3 is related to wind 
utilization for power production. 

With regard to wind utilization, the maximum power that can be extracted by a wind- 
mill depends on the cube of the wind velocity. Hence, in a location where there tend to 
be large steady winds, it would be of advantage to place a windmill such that the ground 
terrain will provide additional localized velocity increases. Possible locations would be 
at the top of a cliff or ridge facing the wind, the top of a circular or oval hill, or in the 
notch between two hills. Except for flow across a two-dimensional cliff or ridge, the 
prediction of the wind behavior involves three-dimensional calculations which are gen- 
erally quite difficult. 

A simplification that will allow some analytical results to be obtained is to consider 



the flow to be inviscid. According to Sutton (ref. 4) the main features of the wind over 
the portion of a hill facing into the wind can be found from an inviscid solution. The in- 
viscid results should provide a good indication of the relative effect on the flow of vari- 
ous types of terrain features. A small perturbation technique will be used to solve the 
inviscid relations, and this provides a limit on how steep the hills can be for the anal- 
ysis to remain valid. For a two-dimensional ridge, all the flow must pass over the 
crest of the ridge. For a hill, however, some of the flow passes around the sides and 
thus provides what is called the three-dimensional relief effect. This effect reduces the 
maximum flow acceleration, and for the same terrain cross-sectional profile it causes 
the perturbation method to be more accurate for the three-dimensional case than for two 
dimensions (ref. 5). 

In the perturbation theory the flow is determined by solving Laplace's equation for 
the velocity perturbation potential. Since this is a linear equation, superposition tech- 
niques can be used to build up solutions. Some solutions will be obtained by using the 
method employed by Scorer (refs. 6 and 7) wherein the flow is first obtained for a two- 
dimensional ridge. Then, results for oval mounds and circular hills can be obtained by 
superposing the ridges at various angles to the incident flow. Some solutions are also 
obtained by using Green’s functions to solve the Laplace equation for the perturbation 
potential. The use of superposition also permits determining the effect of two adjacent 
hills as long as the second hill is not in a region where flow separation would be induced 
by the first hill thereby making the inviscid solution invalid. 


ANALYSIS 


Small Perturbation Equations 


The flow across irregular terrain is idealized by assuming it to be incompressible 
and inviscid. The undisturbed flow is uniform and hence is irrotational; then, as a 
consequence of the inviscid assumption, the flow field remains irrotational everywhere. 
For irrotational flow the velocity components obey the relations 


9v _ 3u _ 0 

dx 3y 


3w 

3y 



3u _ 3w _ q I 
3z 3x J 


( 1 ) 
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(All symbols are defined in appendix A. ) These conditions are satisfied if there is a 
potential function # such that 
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9x 
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3z 




For steady incompressible flow the continuity equation is 

du 9v + 9w _ q 
9x 9y dz 

Inserting equation (2) in equation (3) yields 

9 2 $ t 9 2 <fr | 9 2 # _ Q 
9x 2 9y 2 9z 2 


so that the potential function is obtained by solving Laplace’s equation. 

The local variations in terrain height are assumed sufficiently small so that only 
small perturbations are produced on the incident uniform velocity. Then 

u = Uoo + u’ ^ 

v = v' >■ 

w = w’ > 


Let the velocity potential be 


$(x, y, z) = u TO x + <p(x, y, z) 


where <p(x, y, z) is the perturbation potential such that 



i£ = u * 

3x 

i£ = V’ 

3y 




dcp _ 
3z 


w 




Substituting equation (6) into equation (4) shows that the perturbation potential also satis- 
fies the Laplace equation 


~2 .2 2 2 

i_£ +i_£+ Lse. = o 

, 2 -,2 -.2 

3x 9y 0z 


( 8 ) 


At the surface of the ground the flow streamlines are tangent to the surface; there- 
fore, at the ground 


U oo + u> 

dx/dt 


V* _ w' 
dy/dt dz/dt 


If the approximations of small perturbation theory (ref. 5) are used. 


w l = u 


oo 


dz 

dx 


( 9 ) 


at the surface of the ground (z = h) which gives the boundary condition for <p as 

= Uoo (10) 

0Z 0X 


For small perturbation theory, a Taylor expansion of quantities away from the base 
plane z = 0 shows that to a first-order approximation the boundary condition can be ap- 
plied at the base plane z = 0. Hence, the boundary condition for equation (8) becomes 


i£ = u 

0Z 2=0 00 0x 


( 11 ) 


where h is the local terrain height. 
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Superposition of Ridges to Obtain Velocities for Circular and Oval Hills 


Flow normal to a two -dimensional ridge . - For flow across a two-dimensional ridge 
as shown in figure 1 a solution for (p can be tried in the form of a product solution 

<P = u^fCxJe' 152 (12) 


where k>0 as the perturbation would not increase with height above the ground. Sub- 
stituting equation (12) into equation (8) gives the equation for f(x), 



+ k 2 f = 0 


(13) 


which has the general solution 

f = A sin kx + B cos kx 


Then, from equation (12) 


(p = u co (A sin kx + B cos kx)e~ kz 

Since this is valid for any k > 0, the solutions can be superposed to give 

/*00 

<p = Uoq / [A(k)sin kx + B(k)cos kx] e -kz dk 
* / k=0 


(14) 


Then by differentiating 


u’ = -^ = u f k[A(k)cos kx - B(k)sin kx] e~ kz dk 
9x • / k=0 


(15a) 


w 


’ = = U f (-k)[A(k)sin kx + B(k)cos kx] e" kz dk (15b) 

dz *^=0 


To obtain the equations of the streamlines, equations (15) are substituted into equation (9) 
to yield 
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dz 


/ oo 

k[A(k)cos kx - B(k)sin kx] e' 

:0 


•kz 


dk> = 


/ OO 

(-k)[A(k)sin kx + B(k)cos kx] e~ kz dk 

:0 


This can be integrated to obtain along any streamline 


= f (A(k)cos kx - B(k)sin kx] e dk + z 
*^k=0 ° 


(16) 


where z Q has a different constant value along each streamline. 

To evaluate the coefficients A(k) and B(k) the small perturbation approximation 
that the boundary condition at the surface of the ground can be applied at the base plane 
z = 0 is used. The perturbation theory requires that the ground contour be smooth so 
that the effect of higher harmonics (large k) is small. Then if z Q = 0 is allowed to 
correspond to the ground streamline, equation (16) yields 


h(x) = 



cos kx - B(k)sin kx] dk 


For a specified h(x) the A(k) and B(k) can be found by expanding h(x) in a Fourier in- 
tegral (ref. 8): 


A(k) 


=- r h(D 

TT •'£=-«’ 


cos k| d| 


B(k) = - - f h(£)sin k£ d£ 

jt J%=- oo J 


(17) 


If the same approximation of referring the ground surface to the base plane is used, 
equations (15) yield for a distance 6 above the ground 


u = u. 



k[A(k)cos kx - B(k)sin kx]e _ls ^ dk 


(18a) 
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(18b) 


w = -u 


J ~ k[A(k)sin kx + B(k)cos kxTje -k ^ dk 


Now consider a characteristic ridge contour that has been used in the literature 
(refs. 6 and 7): 


Mx) _ 


h 


m 



(19) 


where b is the value of x at the half-height of the ridge (fig. 1). The h(x) is an even 
function, so from equation (17) B(k) = 0, and the integration for A(k) yields 


A(k) - h m b e" kb 


( 20 ) 


Then from equation (18a) the u becomes 


u 



1 + h m b f k e"^ b+6 ^ k 

m *'k=0 


cos kx dk 


which integrates to 



At the peak of the ridge x = 0, and at height 6 above the peak the velocity is 



( 21 ) 


( 22 ) 


Undisturbed flow at an angle to two-dimensional ridge . - In figure 1 the centerline of 
the ridge was along the y-axis; now let it be at an angle y to the y-axis as shown in 
figure 2. If equation (19) is used, the height of the ridge at location (r, 0) is 


7 



I Hill I 
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h(r, 0) _ 1 1 

2 

1 + ^ cos 2 (0 + y) 

The undisturbed velocity u^ has a component u M sin y parallel to the ridge that is un- 
changed by the ridge. The component u^ cos y normal to the ridge is influenced by the 
ridge as given by equation (21). Then at (r, 0), adding the components in the x-direction 
gives 






which simplifies to 



1 + h b cos 2 yi — < b + 6 > 2 - r2 cos2 (g ^ 


.2 2 2 / 

) + r cos ( 


(24) 


Flow across an oval mound. - Consider an oval mound as in figure 3 with its major 
axis at an angle u> to the y-axis as shown. As discussed in reference 7 a superposi- 
tion of solutions for two-dimensional ridges at various y angles can be used to build up 
the solution for an oval mound. The oval shape is obtained by introducing an angular 
weighting factor <3?(y, cu) into the superposition. The height of the mound is found by 
using equation (23) as 





I 



A typical weighting function is 


#(y, co) = cos 



(26) 


where 


a = y - u> - tt - < (y - co) 

2 

a = y - (jo - — < (y - co) ^ — 

2 2 

Q' = y- co+7T (y - co) < - — 

2 

Although the velocity distribution is a function of angle co, the contour lines of constant 
ground elevation do not depend on co. Hence, equation (25) is most conveniently eval- 
uated for co = 0. This yields 



2 


TT 

Results for the mound height were evaluated by numerical integration for various r/b 
and 9 and will be given later. 

The quantity b in equation (27) is the half-width at the half-height of all the compo- 
nent ridges used in the superposition to build up the oval mound contour. The b must 
now be related to a characteristic dimension of the mound. The characteristic dimen- 
sion will be taken as r Q , which is in the cross section A-A of the minor axis as shown 
in figure 3. This would be a known quantity from the hill geometry. Since co = 0 in 
equation (27), the locations along section A-A are along the x-axis corresponding to 
0 = 0. Then, with h = h m /2 and r = r Q , equation (27) gives 


9 




( 28 ) 


The particular value for r^b that satisfies this equation will be given when numerical 
results are presented. This numerical value fixes the b for a r Q given by the specific 
hill geometry, and hence is used to relate the results in equation (27) to the physical 
geometry. 

To obtain the flow across a mound having the orientation shown in figure 3, the ve- 
locity distribution in equation (24) is superposed by using the same superposition proce- 
dure as was used in equation (25). The result is 



After subst "''ting equation (26) and integrating the denominator, equation (29) becomes 



At the top of the mound r = 0 and equation (30) becomes 
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= 1 + 


7r 2 h b 
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3t/2 

6(u 2 - 8) (b + 6) 2 2) 


r 

•/v=“ 


&(y, U>)cOS y dry 



( 31 ) 


The integral in equation (31) was evaluated numerically for the <f> given by equation (26): 
some results for the peak velocity will be given later as a function of the angle of inci- 
dence u> of the flow relative to the mound major axis. 

Flow across a circular hill. - For the special case of a circular hill (figs. 4 and 5), 
the weighting function 4> = 1 and equation (25) becomes (note h is now only a function 
of r) 



This is integrated to give 

h(r) _ 1 

v "iRj‘ 

At h/h m =1/2, r s r Q so 1/2 = 1 + (r Q /b) 2 which gives b = r Q /y[z. Then the con- 

tour of the circular hill is expressed in terms of r , which would be known from the hill 
cross section, as follows: 



From equation (29) if $ = 1 and r = 0 there is obtained at the top of the axisym- 
metric hill 
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i , h m b 1 f T 

u m * (b + 6) 2 *^=- 


ff / 2 2 ■. 

COS y ay 

r=-(?r/2) 


Integrating and substituting b = r give 


— = 1 + 


V3h 


m 


u. 


m 


2r 


[i + Vi A 

r„; 


( 33 ) 
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Analysis for Circular (Axisymmetric) Hills by Use of Green’s Functions 

General relations . - The perturbation potential cp is governed by the Laplace equa- 
tion (eq. (8)) in the upper half space subject to the boundary condition that the normal 
derivative of <p be specified at the boundary z = 0 as given by equation (11). This is 
known as the Neumann boundary condition, and the solution for cp can be found from 
Green’s fundamental solution as (ref. 9) 

cp(x, y,z) = - — 

Ait 

where g is the Green's function for the 
serted, equation (34) becomes 


// 

x-y plane 


g d| d ri 
3? 


(34) 


Neumann problem. If g (from ref. 10) is in- 



Then u' is found by differentiating with respect to x: 


(35) 
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A form, that is a little more convenient for analytical integration is found by integrating 
by parts to yield 


u 


u’(x, y, z) 

2tt 



rjj— —oo ^ — o 


[(* - x) 2 + ( V - y) 2 + z 2 ] ^ 


^ d| d v 

1/2 .,2 


(36) 


Of specific interest here is the u-velocity at a distance 6 above the peak of the hill (at 
x = y = 0): 


u’(0, 0, 6) 



(37) 


Anal y tical relations for three hill contours . - Solutions will now be found for three 
different hill contours. To illustrate the Green’s function technique in comparison with 
the previous analysis, the hill contour that was obtained by superposition of ridges will 
be considered (fig. 5): 



This is differentiated twice with respect to x and then substituted into equation (37). 
The results are put into polar coordinates to yield 



Integrating over 9 gives 
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R dR 


u*(0, 0, A) _ 3h m 
u^, 2hy, 


f 


0 (H 2 + A 2 ) l/2 [(l + 3 E 2p ( 1 + 3R 2 ) 


2jt 


97tR 


3/2 


5/2 


o 

If 1 + 3R is the square of a new variable of integration this can be integrated; the 
u = u M + u' then becomes 


ju_ _ i + Vi 


u„ 



This is in agreement with equation (33). 

A second circular hill contour that will be considered is 



(39) 


(40) 


This is of interest because it has the same profile that was used for the two-dimensional 
ridge (eq. (19)), and the results using equation (40) will provide a comparison between 
the flow acceleration effects by a two-dimensional ridge and a circular hill having the 
same contour. 

Equation (40) is differentiated twice with respect to x, and the results are substi- 
tuted into equation (37). After putting the results into polar coordinates the velocity 
perturbation above the peak of the hill is given by 



Integrating over 6 gives 
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u T (0, 0, A) _ n m 


u 



A 2 ) 1 ^ 2 j^(l + R 2 )^ (l+R 2 )' 


2 jt 


4jtR 


R dR 


2 2 

If R + A equals the square of a new variable of integration, the previous equation can 
be integrated to yield 


2h 


u(0, 0, 8) _ ! | m 


u 


5 0 3=2 

r Q + 2r 0 6 


4 (4 - 62 ) 5/2 J 


cos 


- 1/6 



o/j 


36r: 


(4 - * 2 y 


(41) 


As a third example to determine the effect of hill contour, the following exponential 
shape will be examined (see fig. 5): 


(42) 

where /3 2 = In 2 = 0. 693. Substituting into equation (37) yields, after placing the results 
in dimensionless polar coordinates, 




- (e - ^ 2 ^ 2 - 2/3 2 R 2 cos 2 0 6"£ 2r2 )r dR dS 

2 . A 2 s ) */ 2 7 


(R + A 


) 


In a similar fashion to what was done for the previous two contours, this is integrated to 
yield 


5) = x + 2 p {- ± g ± + ^ 


1 3 


ir 


r o t 2 r o 2 


2 \r„ 


2n 


1 - erfj 


(43) 


Some results from these relations will be given later. 
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Superposition of the Effect of Two Hills 


In figure 6 two hills are shown, and a quantity of interest is the velocity augmenta- 
tion at the notch between the hills. The flow direction is normal to the line between the 
hill centers. Since the Laplace equation for the perturbation potential is linear, the per- 
turbation potential for two hills can be found by adding the potentials for each of the in- 
dividual hills: 


u = — (u B x + <p+ + <p 2 ) = u oo + u i + u o ( 44 ) ^ 

ax |( 

Within the assumptions of small perturbation theory the surface of the combination of 
two hills is taken to be at the base plane (z = 0) as was the case for a single hill. Then 
to obtain the perturbation at the location of interest (point N) in figure 6, the u’ is ob- 
tained at z = 0 for a single hill at a radius m/2 in the plane normal to the direction of 
the flow. From equation (29) with r = m/2, 9 = n/2, and 6=0, this is 



Using cos [(ff/2) + y\ = -sin y and adding the effect of the two perturbations give at 
point N 



For a single hill for which equation (29) is valid there was obtained b = r Q /\^3 so that 
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( 45 ) 



The perturbation velocity can then be written as 


u' 


u 


OO 




(46) 


The quantity I(m/r Q ) was calculated numerically and results for (u'/u OT )(r o /li ) will be 
given later. 

As a second example to show the effect of superposing two hills, results will be 
calculated for two adjacent exponential hills as described by equation (42). At the ground, 
z = 0, equation (36) gives the velocity perturbation as 



a 2 h 


[(?-x) 2 + (T7-y) 2 ] 1/2 ^ 2 


d£ d rj 


From equation (42) 



so that 



= -2h. 


m 



p/r^+y 2 ) 



Let h = /3h/r 0 , | 


/3£/r 0 , etc. ; the velocity perturbation then becomes 
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/ >2 a 2 \ 
_ +»? ) 


[(I - x) 2 + (77 - y) 2 ] 


1/2 


(i - 2i 2 )di ^ 


As shown in figure 6, the locations of interest for use in the superposition of two hills 
are along the y-axis (x = 0) which gives 



Transforming to polar coordinates gives 



1/2 


(1 - 2R 2 cos 2 0)e" R R dR d0 


(r 2 - 2yR sin 9 + y 2 ) 

If the double integral is given by nG(y), this is written as 

— — = j3G(y) 


(47) 


u °° h m 


(48) 


The analytical evaluation of G(y) is given in appendix B. 

To superpose the effect of two hills the location of interest is at 
y= m/2 [y = (/3/2)(m/r o )]. Then at the notch between the two hills equation (44) yields 

— = 1 + 2/3 — g(& — ^ 


u 


r o \ 2r o/ 


An expression for G along with nmnerical values will be given in appendix B, and the 
velocity perturbation will be plotted in the form 
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— — = 2/3G 

u oo h m 


RESULTS AND DISCUSSION 


Flow Across a Ridge 


The contour of the ridge as given by equation (19) is plotted in figure 7 for a few dif- 
ferent values of h m /b, the ratio of maximum height to the half -width of the ridge at 
half-height. From equation (22) the surface velocity at the peak of the hill is 


= 1 + ■ 


so that the most gradual of the three ridges in figure 7 will give u/u ot = 1. 33 at the 
peak. For the most steep ridge shown, u/u^ = 2 at the peak of the ridge, and since the 
velocity perturbation is not small the accuracy is questionable. The accuracy is prob- 
ably satisfactory for the two other ridges shown with h m /b = 1/2 and 1/3. 


Flow Across a Circular Hill 


The contours of three different circular hills are shown in figure 5 for h m / r 0 = 0. 5. 
From equations (39), (41), and (43) the surface velocity at the peak of these hills is 


-^= 1 + C 


where the coefficient C is given in the following table for various hill contours: 
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The largest C, which provides the largest velocity increase, corresponds to the 
profile with the steepest contour near the peak. For the contours shown in figure 5 
where h /r =0. 5, the velocities at the peak are 1. 37 to 1. 43 times the undisturbed 
free stream velocity. Since these perturbations are not small compared to 1 the ques- 
tion arises as to whether results from small perturbation theory are meaningful for 
hills with h m /r Q this large. To obtain an indication of the error, a rough comparison 
can be made with results for flow over ellipsoids (ref. 11). These results indicate that 
the small perturbation theory is in error only about 5 percent for h/r c = 1/3, so the 
error for h/r Q = 1/2 is probably acceptable. For h m /r 0 < 1/4 the theory should be 
quite accurate. 

The second contour in the table of C values is the same as the contour of the two- 
dimensional ridge considered previously. Comparing with equation (22) shows that at 
the peak the acceleration by the hill is 78. 5 percent of that for the ridge. This difference 
is the ’’three-dimensional relief effect" and results from part of the flow going around 
the hill rather than all passing over the crest of the ridge. 


Flow Across an Oval Mound 

The local height of an oval mound is given by equation (27) as a function of a dimen- 
sionless radial coordinate r/b and angle 0 as shown in the inset of figure 8. Contours 
were obtained by numerical integration of equation (27) and are shown in figure 8 for 
various angles where 0 = 0 corresponds to the cross section across the most narrow 
part of the mound. From equation (28) the characteristic length r Q is defined as the 
radius within the cross section 0 = 0 at which the mound height is one-half the maxi- 
mum height. From figure 8 there is obtained (at 0 = 0 and h/h m = 1/2) r Q /b = 1. 22 so 
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that b = 0. 820 r . From the size of the actual mound, the r Q is obtained at the half- 
height of the most narrow cross section; thus, the value of b is then known. For con- 
venience a second abscissa scale of r/r Q is given in figure 8. Figure 9 shows the 
shapes of contours of equal height. The mound is approximately twice as long as it is 
wide. 

From equation (31) the surface velocity at the peak of the mound is 

22 » 2 r n cos 2 yd, 

6(tt 2 - 8 ) 


U -f U. H 

— = 1 + — = 1 + 


u 1 


m 


u 


u 


The quantity in square brackets, which is equal to (r 0 /h m )(u'/u 00 ), was evaluated nu- 
merically and is plotted in figure 10 as a function of the angle co which specifies the 
wind direction. When = 0 the wind is blowing along the short axis of the mound and 
this provides the maximum velocity enhancement: 


— (w = 0) = 0. 905 — 


u 


As the wind shifts away from the direction along the short axis, the perturbation de- 
creases and reaches a minimum for u> = 90° where 


— (w = 90°) = 0. 315 — 

U °o r O 

The contour for 6 = 0 in figure 8 is very close to that for the ridge in figure 7 and 
the l/[l + (r/r Q )^] hill contour in figure 5. The velocity perturbation at the peak of the 
ridge or hill is in all cases equal to a constant times the maximum height of the ground 
divided by the horizontal distance to the location of the ridge or hill half-height. The 
constant is 1 for the ridge, 0. 905 for flow parallel to the shorter axis of the oval 
mound, and 0. 785 for the circular hill. 


Flow In Notch Between Two Adjacent Hills 

A possible location for wind velocity enhancement is at the notch between two hills. 
This geometry is illustrated in figure 6 and the location of interest is at point N. The 
wind direction is normal to the line connecting the centers of the two hills. The coor- 
dinate l is along this line and extends from the symmetry plane of the combined geom- 
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etry of the two hills. The cross section A-A in figure 6(b) is given in figures 11(b) 
and (c) and 12(b) for two hill contours having various spacings apart. 

Figure 11 deals with the superposition of hills having the contour 


m 


ft; 


Part (a) of the figure shows the hill contour for three values of h m / r 0 - As shown by 
figure 6, m is the spacing between the hills. Adding the heights of the hills for two dif- 
ferent spacings yields parts (b) and (c) of figure 11. These two parts show the cross 
section A-A as shown in figure 6. Figure 12(a) shows the contours of two single hills 
having the shape 



The cross section when two of these hills are superposed is shown in figure 12(b). 

For the hill contour in figure 11 the surface velocity perturbation at point N in the 
notch between the hills is given by ^3 I where I is the quantity in square brackets in 
equation (45). This has been evaluated numerically and is plotted in figure 13 as a func- 
tion of spacing between the hills. For the hills shown in figure 12 the velocity perturba- 
tion at the ground midway between the hills is given by equation (49). This has been 
evaluated in appendix B and is plotted as the dashed line in figure 13. 

Figure 13 shows how the velocity perturbation at the notch between the hills dimin- 
ishes as the spacing between the centers of the two hills is increased. In the vicinity of 
r/r 0 = 2 the hill height in figure 12 decreases much more rapidly than for the contour 
in figure 11. As m/r Q is increased in figure 13, this accounts for the rapid drop in the 
dashed curve as compared with the solid curve. The results show that for a typical case 
of hills with h m /r Q = 1 spaced at m/r Q = 4 the velocity perturbation is about 20 percent 
of the undisturbed velocity. For the same spacing, if the hill height is reduced in half so 
that h m /r G = 1/2, then the perturbation is reduced to about 10 percent. 


CONCLUSIONS 


Potential flow solutions have been obtained for wind passing over some forms of 
ground topography that can provide a localized velocity increase. The geometries studied 
are axisymmetric hills, an oval hill, a two-dimensional ridge, and the notch between two 
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circular hills. These are possible terrain features to provide locally increased veloc- 
ities for windmill siting. 

The hills and the ridge have cross sections typical of those found naturally. The 
characteristic hill or ridge dimensions are the maximum height h m and the horizontal 
distance from the center line or center plane to the location of the half-height; this hori- 
zontal distance is r Q for a hill and b for a ridge. A quantity of interest is the velocity 
perturbation u* at the ground at the peak of the hill or ridge where the maximum veloc- 
ity increase occurs from the uniform incident velocity. The velocity perturbation was 
found to equal h m /b for the ridge shape that was analyzed, and it was equal to about 
0. 8 h m /r Q for a hill. The coefficient of 0. 8 for the hill is an average value, as it de- 
pends somewhat on the specific hill shape. A range of 0. 74 to 0. 87 was found for some 
typical hill contours. Comparing an axisymmetric hill with a ridge having the same con- 
tour showed that the increase of velocity at the peak of the hill was 78. 5 percent of that 
for the ridge. This is the three-dimensional relief effect wherein some of the flow 
passes around the hill rather than all passing over the crest as is the case for a two- 
dimensional ridge. 

A superposition of ridges was used to obtain the flow over an elongated mound with 
oval shaped contours of constant height. Calculations were made showing the effect of 
the flow approaching the mound from various directions relative to the mound major 
axis. When the flow was along the short dimension of the mound the behavior was close 
to that across a two-dimensional ridge. A much smaller effect was produced when the 
flow was along the long axis of the mound. This would be expected as the flow is unaf- 
fected when moving parallel to a two-dimensional ridge. 

Superposition was used to obtain the velocity increase at the notch between two hills. 
Increases of 20 percent are possible for hill contours of reasonable shape. Larger in- 
creases can be obtained with hills having greater slopes, but these cannot be calculated 
with the small perturbation theory used here. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 1, 1976, 

506-24. 
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APPENDIX A 
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SYMBOLS 

coefficients that are functions of k 
value of x at half -height of a ridge 
integral defined by eqs. (47) and (48) 

Green's function for Neumann problem 
local height 
maximum height 

integral defined by eqs. (45) and (46) 
number > 0 

coordinate along cross section of two hills (fig. 6) 
spacing between two hills (fig. 6) 
dimensionless coordinate, r/r Q 
radial coordinate 

radial location at which hill is at half-height 
time 

velocity components in x-, y-, z-directions 

velocity component perturbations 

undisturbed velocity approaching irregular terrain 

rectangular coordinates 

dummy variable in eq. (26) 

constant equal to In 2 = 0. 693 

angle of ridge line (fig. 2) 

height above peak of hill 

dimensionless coordinate 5/r Q 

angular coordinate (see figs. 2 or 3) 

coordinate normal to ridge line (fig. 2) 

dummy variables in the x-, y-, and z-directions 

velocity potential function 
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c p perturbation velocity potential function 

u) angle of major axis of mound (see fig. 3) 
Superscript: 

~ variable multiplied by /3/r Q 



APPENDIX B 


EVALUATION OF INTEGRAL G(y) a 
From equation (47) the integral G(y) is 



(1 - 2R 2 cos 2 e) e~^ _ RdRde 


(r 2 - 2yR sin 9 + y 2 ) 


1/2 


In rectangular coordinates, as given by the relation preceding equation (47), this is equal 
to 



( 1 - 2 gV(jV) 
[? 2 + (5 - y ) 2 ] 


d| drj 


Let a new variable be y = t] - y to give 



(1 - 2£ 2 )e^ 2+ fr + y )2] 


(l 2 + r 2 ) 


1/2 


d y 


2 A n 9 A 

Changing to polar coordinates where r = + y and <p = tan (y/|) yields 

*2jt 

. A i 

G(y) = - 



(l-a^eosW^^ 8111 ^ . dr a* 


n a o a 

/, o 2 2 \ -r -y -2ry sin cp 

(1 - 2r cos <p)e J e J ^ 


dr d cp 


a This integration was carried out by William F. Ford of NASA, Lewis Research 
Center. 
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Since G(y) = G(-y), because the sin jp is integrated over an equal range of positive 
and negative values, the factor e -2r y s * n ^ can be replaced by (l/2)(e^ r ^ s * n ^ 

+ e"^ r y sin c f > ) = cosh(2ry sin cp) = cos(2iry sin cp) so that G(y) becomes 


G(y) = -i 


IT 



2 . -r 2 -y 2 
cos <p)e J 


cos(2iry sin cp)Ar d (p 


(Bl) 


From the definite integral tables (ref. 12, p. 143), 


f 

J r=0 


e -r cos(2iry sin <p)dr = e^ s * n ^ 


f 


2 ^2 2 

2r 2 cos 2 <p e _r cos(2iry sin <p)dr = cos 2 <p e^ s * n ^ (2y 2 sin 2 <p + 1) 


The G(y) in equation (Bl) is then 


G(y) = 


/»2 V 

/ 

J(p=0 


i y cos ^ sin 2 <p(l - 2y 2 cos 2 <p)d<p 


2 2 

The double angle formulas are used that cos 9 = (1/2)(1 + cos fi) and sin 9 
(1/2)(1 - cos j3) where /3 = 20 and the G(y) reduces to 


G(y) = 


(i/2)> r 71 r 

^ U 


1 -2 

- cos (3 — y (1 - cos 2/3)1 e 
2 


_ (l/2)y 2 cos /3 


d/3 (B2) 


From reference 12, page 190, equation (B2) is integrated to give 


G(y) = 



[d - y 0 )^ J 0 ^y 0 ) - ^Ji(iy Q ) - yc/V 1 ^ 


where y Q = (l/^y 2 . 


Using the following identities in equation (B3) 


(B3) 
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J o^o ) = W 


iJlUYo) = -^(Yq) 


' J 2 (iy o> - W ‘ 7 W 
y n 


gives G(y) in the final form 

(B4) 

Some values of G(y) evaluated from equation (B4) are 


y 

G(y) 

y 

G(y) 

0 . 00 

0 . 886 

1.40 

0.234 

. 20 

. 860 

1.60 

. 166 

.40 

.787 

1.80 

. 117 

.60 

.679 

2.00 

.0826 

. 80 

. 556 

2.45 

.0409 

1.00 

.433 

2.83 

.0204 

1.20 

.324 

3.16 

.0173 
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Figure 1. - Two-dimensional ridge with local height h(x) in wind with 
uniform velocity u M . 



Figure 2. - Coordinates looking down on two-dimensional 
ridge at an angle to the wind. 



Figure 3. - Coordinates viewing oval mound from above, and 
geometry of cross section along minor axis of mound. 
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Figure 4. - Circular hill having local height h(x,y) 
in wind with uniform incident velocity, Uoo. 




Figure 5. - Geometry for circular hill. 
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^Location of interest, N 
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(b) Coordinates looking down on hills. 

Figure 6. - Geometry for flow across two superposed hills. 
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Height contours of oval mound along radial directions at various angles from the minor axis. 




(a) Contours of single hills. 



llh m 

(b) Superposition of two hills, ml r 0 = 4. 





(b) Superposition of two hills, ml r 0 = 4. 

Figure 12. - Superposition of two hills each having contour h/h m = exp [-(In 2)(r/r 0 ) 2 ], 
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